Nuclear spin relaxation rates in two-leg spin ladders 
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Using the transfer-matrix DMRG method, we study the nuclear spin relaxation rate 1/Ti in the 
two-leg s=l/2 ladder as function of the inter-chain (J±) and intra-chain ( Jy ) couplings. In particular, 
we separate the q y = and n contributions and show that the later contribute significantly to the 
copper relaxation rate 63 (1/Ti) in the experimentally relevant coupling and temperature range. We 
compare our results to both theoretical predictions and experimental measures on ladder materials. 
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Dynamical properties, while providing the most de- 
tailed information on the physics of low-dimensional an- 
tiferromagnets are also the most difficult for theoretical 
predictions. Measured in NMR experiments probing the 
low-frequency spin dynamics, they continue to reveal un- 
expected features. For instance, the recent experiments 
by Imai et al. (!]] on the CU2O3 two-leg s= 1/2 ladders 
in La6CagCu2404i have shown the different tempera- 
ture T dependence of the nuclear spin relaxation rate 
1/Ti for the oxygen ( 17 (l/Ti)) and copper ( 63 (1/Ti)) 
atoms. Although both exhibit activated behavior below 
T ~ 350 K, 63 (1/Ti) shows a crossover to a linear T de- 
pendence above T ~ 350 K. Shortly after these experi- 
ments were performed, Ivanov and Lee [Q proposed that 
the crossover may be understood from the momentum 
transfer q = (jr, it) contributions in the weakly coupled 
chain limit. The importance of the q y = tt processes was 
also mentioned earlier by Sandvik et al. Q in a study of 
63 (1/Ti) for SrCu2C>3, where they consider the isotropic 
ladder at T/J > 0.2. Despite these findings, a better 
theoretical understanding of 1/Ti in isolated ladders is 
necessary, as emphasized in Refs. |Qp). 

Our purpose is to determine 1/Ti for the experimen- 
tally relevant ratio of inter (J±) and intra-chain (Jy) cou- 
plings, separating the q y = and tt contributions. We 
will use the transfer-matrix density-matrix renormaliza- 
tion group for the evaluation of thermodynamic proper- 
ties j| and local imaginary time correlations ||||] of one- 
dimensional tight-binding models. Combined with an an- 
alytical continuation, this technique has recently proven 
to provide reliable real frequency correlations down to 
low temperatures In the low-T regime, our results 
show a good agreement to analytical results obtained in 
the weak or strong coupling limits [||j7]]. At higher T, we 
demonstrate the importance of the (1/Ti) ff contributions 
to 63 (1/Ti) experiments, in particular, our result repro- 
duces the crossover to the linear paramagnetic regime 
observed in the 63 (1/Ti) rate. 

The most recent estimates of J\\/J± from susceptibil- 
ity and neutron scattering experiments agree for a value 
of about 0.5 (HJ9). On the other hand, estimating J\\/Jj_ 
from NMR rates, we found that the standard ladder, to- 
gether with the hyperfine couplings measured in Refs. 



PIPI, are unable to account for J\\/Jj_ ~ 0.5 but favors 
J\\/J±~l. 

The two-leg ladder is described by the antiferromag- 
netic Heisenberg Hamiltonian 
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where S nt i denotes a s=l/2 spin operator at the i-th rung 
and the n-th chain. It is by now well established that the 
spectrum of the Hamiltonian ([!]) consists of an S = 
ground state, the lowest-lying excited states forming a 
gaped S = 1, k y = tt single magnon branch with mini- 
mum A at k x = tt J^lOpr]] . 
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FIG. 1. Comparison of the susceptibility and Knight shifts 
to transfer-matrix DMRG results for J\\/J± = 0.5 and 1. In- 
set: susceptibility and Knight shift measurements compared. 

The issue of determining the values of J11 and Jj_ has 
raised some controversy, indeed estimates from various 
authors jl 12 ,13 MO range from J±/J\\ — 0.5 to 1, with 
an emerging consensus for about 0.5. Fitting our numer- 
ical results to susceptibility ]l3] and Knight shift [Q ex- 
periments, we confirm in FigTll^ that the best agreement 
is found when J±/J\\ = 0.5 for SrCu203 15 1, implying 
energy scales A ~ 420K and Jy ~ 2000K. Nevertheless, 
a reasonable fit to J±/J\\ = 1 is also possible in a slightly 
narrower T— range, showing that the susceptibility is rel- 
atively insensitive to J±/J\\ at low— T. Moreover, we 
point out in the inset that the susceptibility and Knight 
shift measurements are not fully consistent through the 
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entire T— range, implying material-dependent Jm values. 
In principle, Knight shifts are better suited for a compar- 
ison with theoretical results, since they involve no sub- 
traction of a Curie term and no unknown g factor |l2||§ . 
The nuclear relaxation rate 1/Ti is given by 
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J b F(q)S(q,u; N ) dq x , 6 = 63,17 (2) 

q y =0..TT 

S ^^ = Jn ^ J Vn, i (*)Sm, J )e i(q " W - i,m - n)+ut) d* • 

Here, con is the nuclear Larmor frequency and b F(q) are 
the appropriate hyperfine couplings. According to Ref. 
@, 63 F(q) = A 2 for 63 Cu, and 17 F 2 (q) = A F 2 cos 2 (q y /2) 
for the rung 0(2) oxygen atoms |ll|. Contrary to CuC>2 
planes, it has been argued that considering only a lo- 
cal hyperfine interaction for the 63 Cu nucleus in ladders 
P,[l3"| is sufficient. 

To determine 1/Ti, the transfer-matrix DMRG 
method permits a very precise evaluation of the imag- 
inary time Green's function j(J 



G(r) = (Sf^S^j) 
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where S*At) 
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and the indices (m,j) £ 



{(1, 1), (2, 1), (1, 2), (2, 2)} run over the four corners of a 
plaquette. For all calculations, we chose an imaginary 
time slice e = 0.025/Jy and m = 200 states were kept to 
represent the transfer-matrix. After an analytical con- 
tinuation using the Maximum Entropy method, we can 
resolve S(q,u>) in q y = or tt and obtain the averages 
over the q x momentum transfer: 



S (q y ,uj) = — J dq x cos (q x /2)S(q x ,q y ,uj) 



S-„{q yi uj) = — J dq x sin (q x /2)S(q x ,q y ,uj) . (4) 

As indicated by the notation, So(q y ,w) (S w (q y ,u>)) is 
dominated by processes with q x close to (tt). For the 
hyperfine couplings given above, the copper and oxygen 
rates are fully determined by appropriate combinations 
of Sq^fayiOj), q x — 0, it. Accordingly, we also define 



S (q y ,u N ) + S„{q v ,UJ N ) 



(5) 



the dimcnsionless contributions to 1/Ti from the q y — 
0, 7r momentum space sectors. 

For illustration, we show in Fig. || Sq x (qy,u>) for the 
case J±/J\\=l as a function of T. As the latter drops, we 
can see clearly the signatures of the low-lying spectrum 
emerging. Indeed, the main peak is due to excitations 
from the ground state to the single magnon branch for the 
q y =Tr sector and to the 2— magnon continuum for ^=0, 
for instance, it appears correctly from S^itr, to) that the 



minimum gap approaches the T=0 value A=0.502. Such 
contributions to Sq !c {qy,Lo) become T independent for 
T^O, however, they do not contribute to 1/Ti as they 
all involve frequencies ui>A~^>lun (~ 3 mK). On the 
contrary, there are thermally activated contributions to 
1/Ti, as they involve excitations for cu— >0 |l7[ . As seen 
in Fig. ^, this limit is dominated in the low-T regime by 
processes with (q x , q y )={0, 0) and (tt,tv). 
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FIG. 2. S qx {q y ,uj) for J ± /J ]} =\. T/J r 
1/3 (dashed), 1/4 (thin) and 1/6 (thick). 



1/2 (dotted line), 



The first contributions involve two thermally ex- 
cited magnons near the minimum of the branch 
e(k x , k y =ir)=A+a(k x — tt) 2 , according to Ref. 0, they 
lead to 



1 

Ti 



1 



' A/T [0.809 - ln( 



T ' 



(6) 



Corrections to the quadratic minimum |l8[ only change 
the prefactor of the exponential term slightly for T/A < 
0.3. The logarithmic divergence in Eq.@ cannot be re- 
solved by the Maximum Entropy method, however, it 
does not change the dominant exponential behavior. For 
instance, the factor in the square brackets changes only 
by about 5 percent from T = A/2 to A/4. 

The second important contributions with q y — tt in- 
volve scattering of single magnons with the 2— magnon 
continuum. By representing the low-lying excitations in 
terms of free massive fermions Jll| in the weakly coupled 
chain limit, it was proposed in Ref. that 



-2A/T 
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As the scale 2A corresponds to the bottom of the 2- 
magnon continuum at k = (0,0) independent of Jj_/J\\ 
|l0|,[l9]], this result should hold over a wider coupling 
range. In fact, exact diagonalisation (ED) results show 
the existence of large matrix elements |(n|5||m)| 2 be- 
tween the 2— magnon continuum at k = (0,0) and the 
single magnon branch. Such processes are characterized 
by a momentum transfer q x which rapidly shifts close to 
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7r with decreasing J±. From the analysis of the spec- 
trum 0, it follows that q x /n w 0.5,0.8 and 0.95 for 
Jj-/J r || = 2, 1 and 0.5. Other important q y = n processes 
occur between the single magnon branch at k=(0, tt) and 
the continuum near k=(7r,0) when J±/Jh < 0.5. How- 
ever, these have a larger activation gap of about 3. 7 A. 

Let us first discuss our results for 1/Ti as a func- 
tion of J± and T. We should point out that the values 
Sq x (q y ,uj — > 0) become very small compared to the mag- 
nitude of the main peak as T is lowered. Therefore, for a 
reliable estimate of 1/Ti, it is essential, (i) to work with 
very precise imaginary time data, and (ii) to separate the 
different q space contributions, as the integrated weight 
of the q y = and tt sectors differ considerably. In Fig. 
||, we present (1/Ti) and (l/Tx) OT in the low-T regime, 
where T is scaled by the spin gap A. 
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FIG. 3. Low-T behavior of (1/Ti) and (l/Ti)„ for 
Jx/Jy = 5 (A), 2 (>), 1.5 (+), 1 (o), 0.8 (□) and 0.6 (O). 
Insets: The dot-dashed lines have slopes 1 (2) for q y =0 (tt). 

(1/Ti)o shows the correct behavior of Eq.(|6|), especially 
for large J±/J\\, when low enough T/A can be reached 
to clearly identify the activated regime. This is verified 
in the inset where we have plotted a line of slope one on 
a logarithmic scale. We also observe that the prefactor of 
the exponential strongly depends on J±/J\\, consistently 
with the factor l/a in Eq.(^jJ). Indeed, as the dispersion 
of the single magnon branch flattens with increasing J± , 
l/a grows and liiaj ± ^ oc l/a = oo [[7|JTo|] . The q y = tt 
results reveal a larger gap which can be fitted consis- 
tently to Eq.(0) when J±/J\\ = 1.5,2 (see inset). When 
Jj_/J\\ < 1, we find some deviations from Eq.(^) In this 
regime, the low-T extraction of (l/Ti)^ becomes more 
delicate as the growing main frequency peak near lu = A 
(due to excitations from the ground state to the single 



magnon branch) may bias the limit S^q, ui — ► 0) (see 
Fig. ||). On the other hand, we cannot exclude that 
the above mentioned processes with a gap ~ 3.7 A con- 
tribute significantly at our lowest temperatures. For the 
case J±/J\\ = 5 (not shown in Fig. ||), (l/Ti)^ = since 
the single magnon band and the 2— magnon continuum 
do not yet overlap. 
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FIG. 4. (1/Ti) and (l/Ti) T as function of T/J\\ for 
J]_/J|j = 1.5 (+), 1 (o), 0.8 (□), 0.6 (O) and 0.4 (<]). In- 
set: comparison of S(0, tt, ui) + S(tt, tt, ui) from ED and Trans- 
fer-Matrix DMRG. 

A wider T range is shown in Fig. |[ (I/Tl)^ obtained 
at T/J||= 1 agree within 5 percent with those from ED 
of an 8-rung ladder, as seen in the inset for Jj_/J||= 0.6. 
However, due to finite size effects, the determination of 
1/Ti from ED becomes meaningless below T/J\\ ~ 0.8. 

As (l/Ti)o obeys a predominantly linear intermediate 
to high-T behavior for all J± above T ~ 0.4 Jy, (l/Ii)-* 
exhibits a qualitative change as J±/J\\ is decreased. In- 
deed, for J±/ J\\ < 0.6 (1/Ti)^ develops a maximum near 
Tyj||~ 0.2 becoming sharper with decreasing inter-chain 
coupling. We believe such a feature is reminiscent of 
the weak coupling limit, as our (1/Ti) ff for Jj_/J|| < 0.6 
shows a good agreement to the behavior predicted in Ref. 
[^J. There, Ivanov and Lee also argued that with increas- 
ing T, (1/Ti) w flattens to the single Heisenberg chain 
result p3[ , as we observe for Jj_/J\\ — 0.4 and 0.6 above 
T/J\\ — 0.4. In spite of the larger gap, we find that 
(1/Ti) ff dominates over (l/Ti)o in a wide intermediate- 
T regime as can be verified from Figs. || and ^. For 
instance, when T/A = 1, (l/Ti)^ exceeds (1/Ti) by a 
factor 2.5 (resp. 9) for J±/J\\ = 1 (rcsp. 0.6). In partic- 
ular, for J±/J\\ = 1, we observe the crossing of (l/Ti)^ 
and (l/Ti)o at T/A ~ 0.35. This consideration implies 
that both q y — and tt contributions are relevant to cop- 
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per 63 (1/Ti) experiments above T ~ 0.35 A (typically 
w 200 K in CU2O3 ladders). In fact, taking into account 
the usually omitted q y = it processes to fit the low-T 
63 (1/Ti) data of Ishida et al. (l| on SrCu 2 3 , we obtain 
for T between 100 and 300 K a gap A « 520 K, while 
considering only the form for (l/Ti)o leads to A « 700 K. 
Hence, the discrepancy between estimates of the gap from 
63 (1/Ti) and susceptibility (leading to A ~ 420K) 

is reduced. 
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FIG. 5. Comparison of NMR rates in SrCu 2 3 @ and 
La6CasCu2404i jl[ to numerical results. Inset: Imai's data, 
the high-T regime extrapolates to zero. 

Let us now compare our results with experimental mea- 
surements of 63 (l/7i) and the rung oxygen 17 (l/Ti) rate 
in SrCu203 and La6Ca 8 Cu24 04i. To convert our val- 
ues to the experimental unit, we considered the hyperfine 
couplings given in Ref. [[[3) for the 63 Cu and in Ref. jjj for 
the 17 nuclei, so that the remaining free parameters are 
only Jj_ and J|| , which we have chosen according to those 
obtained from the susceptibility fits in Fig. |l|. Consider- 
ing that T scales proportionally and 1/T± inversely pro- 
portional to Jii , we found good overall magnitudes, how- 
ever, the precise T-dependence is not reproduced very 
accurately. Our results tend to indicate that J±/J\\ ~ 1, 
in disagreement with J±/J\\ ~ 0.5 from the susceptibil- 
ity. Especially, the 'linear' high-T behavior of 63 (1/T!) 
extrapolating to zero (inset) cannot be reproduced with 
a smaller ratio J±/Ji\. 

These observations raise two important issues about 
ladder materials. First, the hyperfine couplings measured 
in Refs. jj];[l3) may not be sufficient for a quantitative 
discussion of NMR experiments within the standard lad- 



der model. A possible improvement may be to consider 
transfered fields 63 T oc (1 + 2rcosq x ) 2 + r' cos 2 (q y /2) 
on the copper atom, which can suppress the (n, n)- 
fluctuations by large factors and favor a smaller ratio 
J±/Jn. This suggests that a better knowledge of the hy- 
perfine couplings in ladder materials is needed for a quan- 
titative interpretation of measured NMR data. Second, 
small corrections to the standard ladder such as interlad- 
der or frustration couplings are known to have little effect 
on bulk properties such as the spin gap (^3|. Therefore, 
the accessible susceptibility measurements (T/A < 1) are 
relatively insensitive to such corrections. On the other 
hand, local dynamic properties as NMR rates can be 
more substantially affected. 

To summarize, we have shown that our technique pro- 
vides reliable NMR rates 1/Ti in the standard ladder over 
a wide temperature range. In particualar, we demon- 
strated that the q y = n contributions are crucial for un- 
derstanding the crossover observed in 63 (1/Ti) experi- 
ments. Finally, we point out how quantitative results 
can indicate inconsistencies in estimates of J±/J\\ from 
suseptibility and NMR measurements. 
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